Expansions for Quantiles and Multivariate Moments 

of Extremes for 

Distributions of Pareto Type 

by 

Christopher S. Withers 
Applied Mathematics Group 
Industrial Research Limited 
Lower Hutt, NEW ZEALAND 

Saralees Nadarajah 
School of Mathematics 
University of Manchester 
Manchester M13 9PL, UK 



Abstract: Let X nr be the rth largest of a random sample of size n from a distribution 
F(x) = 1 - J2Zo dx~ a ~ il3 for a > and > 0. An inversion theorem is proved and used to 
derive an expansion for the quantile and powers of it. From this an expansion in powers 

of (n _1 , n~P/ a ) is given for the multivariate moments of the extremes {X nin _ s ., 1 < i < k] /n l / a 
for fixed s = (si, . . . , Sfc), where k > 1. Examples include the Cauchy, Student t, F, second 
extreme distributions and stable laws of index a < 1. 

AMS 2000 Subject Classification: Primary 62E15; Secondary 62E17. 

Keywords and Phrases: Bell polynomials; Extremes; Inversion theorem; Moments; Pareto; 
Quantiles. 



1 Introduction and Summary 



For 1 < r < n, let X nr be the rth largest of a random sample of size n from a continuous 
distribution F on R, the real numbers. Let / denote the density of F when it exists. The 
study of the asymptotics of the moments of X nr has been of considerable interest. McCord 
(1964) gave a first approximation to the moments of X n \ for three classes. This showed that a 
moment of X n \ can behave like any positive power of n or n\ = log n. (Here log is to the base 
e.) Pickands (1968) explored the conditions under which various moments of {X n \ — b n )/a n 
converge to the corresponding moments of the extreme value distribution. It was proved that 
this is indeed true for all F in the domain of attraction of an extreme value distribution 
provided that the moments are finite for sufficiently large n. For other work, we refer the 
readers to Polfeldt (1970), Ramachandran (1984) and Resnick (1987). 
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The asymptotics of the quantiles of X nr have also been studied. Note that U nr = F(X nr ) 
is the rth order statistics from £7(0, 1). For 1 < T\ < T2 < ■ • • < < n set U n>T = {U nr ., 1 < 
% < k}. By Section 14.2 of Stuart and Ord (1987), U n has the multivariate beta density 

k 

U n , r ~ B(u : r) = 11 (u i+1 - Ul ) r ^-^- 1 /B n (v) (1.1) 
i=0 

on < u\ < ■ • • < Uf. < 1, where no = 0, u^+i = 1, tq = 0, r^ +1 = n + 1 and 

k 

B n (r) = l\B(r i ,r l+1 -r i ). (1.2) 

i=i 

David and Johnson (1954) expanded X nri = F~ l (JJ nTi ) about u n i = EU nri = Ti/{n + 1): 
x nn = ^2jLoG^(u n i)(U n i — u n i)^/j\, where G(u) = F _1 (u), and using the properties of (jl.ip 
showed that if r depends on n in such a ways that r/n ->p£ (0,1) as n — > oo then the mth 
order cumulants of X n>T = {X nn , 1 < i < k} have magnitude 0(n 1_m ) - at least for n < 4, 
so that the distribution of X n ^ r has a multivariate Edgeworth expansion in powers of n -1 / 2 . 
(Alternatively one can use James and Mayne (1962) to derive the cumulants of X n ^ r from those 
of U n>r .) The method requires the derivatives of F at {i 7-1 ^), 1 < i < k} so breaks down if 
Pi = or pk = 1 - which is the situation we study here. For definiteness, we confine ourselves 
to F~ l {u) having a power singularity at 1, say F~ l {u) ~ (1 — u)~ l l a as u — * 1, where a > 
that is, 

1 - F(x) ~ (1.3) 
as x — > oo. For a nonparametric estimate of a see Novak and Utev (1990). 

Distributions satisfying (|1.3p are known as Pareto type distributions. These distributions 
arise in many areas of the sciences, engineering and medicine. Some of these areas - where 
publications involving Pareto type distributions have appeared - are: hydrology, physics, wind 
engineering and industrial aerodynamics, computer science, water resources, insurance math- 
ematics and economics, structural safety, material science, performance evaluation, queueing 
systems, geophysical research, ironmaking and steelmaking, banking and finance, atmospheric 
environment, civil engineering, communications, information processing and management, high 
speed networks, lightwave technology, solar energy engineering, supercomputing, natural haz- 
ards and earth system sciences, ocean engineering, optics communications, reliability engineer- 
ing, signal processing and urban studies. 

In Withers and Nadarajah (2007a) we showed that for fixed r when (|1.3p holds the dis- 
tribution of X njn i_ r (where 1 is the vector of ones in 3f? fc ), suitably normalized tends to a 
certain multivariate extreme value distribution as n — > oo, and so obtained the leading terms 
of the expansions of its moments in inverse powers of n. Here we show how to extend those 
expansions when 

oo 

F-\u)=Y,bi{l-ur (1.4) 
i=0 
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with ao < ol\ < ■ ■ ■ , that is, {1 — F(x)}x^ 1 ^ a ° has a power series in {x~ Si : Si = (a, — a,o)/oto}. 
Hall (1978) considered (|1.4p with oti = i — 1/a, but did not give the corresponding expansion 
for F{x) or expansions in inverse powers of n. He applied it to the Cauchy. In Section 2, we 
demonstrate the method when 

oo 

l-F(x) = x- a J2 c i x ~ i ^ ( L5 ) 

8=0 

where a > and (3 > 0. In this case, (jl.4j) holds with a, = — l)/a. In Section 3, we 
apply it to the Student t, F and second extreme value distribution and to stable laws of 
exponent a < 1. Appendix A gives the inverse theorem needed to pass from (jl.5p to (jl.4p . 
and expansions for powers and logs of series. 

We use the following notation and terminology. Let (x)i = T(x + i)/T{x) and < x >i= 
T(x + 1)/T(x — i + 1). An inequality in consists of k inequalities. For example, for x in 
C fc , where C is the set of complex numbers, i?e(x) < means that Re(xi) < for 1 < i < k. 
Also 1(A) = 1 or for A true or false and 5%j = I(i = j). For G C k let 6 denote the vector 
with 6i = Y^j=i 8j- 

2 Main Results 

For 1 < ri < • • • < r/% < n set Si = n — r{. Here, we show how to obtain expansions in inverse 
powers of n for the moments of the X n ^ s for fixed r when (jl.4p holds, and in particular when 
the upper tail of F satisfies (|1.5D . 

Theorem 2.1 Suppose 1 1.5]) holds with c$, a, f3 > 0. Then F (u) is given by (fi.^[ ) wrai/i 
aj = ia — 1/a, a = f3/a and b{ = C it i/ a , where Ci^ = Cq Q(— ip, cq, a;* o/ (E32P anrf x* = 
iJ(a,l,c)o/ra): 

/-r _ ip 

= ■0c[('" a ~ 1 ci, 

= Vcr 2a ~ 2 {c c 2 + (V-2a-l) C 2/2}, 
C 3V; = V^" 3a " 3 [cqC 2 + (^- 3a- l)c cic 2 + {(^ + 1) 2 /6(V + 3a/2)(a + 1)} c?] , 
and so on. Also for any 9 in 

oo 

{F-\u)} = J2^-^ ia '^ (2.6) 

at ip = 0/a. 

Note 2.1 On those rate occasions where the coefficients di = C it i/ a in F -1 (it) = Yli=o(^ ~ 
u ya-i/a^. are fc nown f rom SO me alternative formula then one can use = d d Q Ci(0, l/do,d) 

of fSM). 
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Proof of Theorem 2.1 By Theorem A.l with k = 1, u = x a , x = c, we have x a 
Ei=o< (! - M ) 1+ia at n = where 



^2 = c o 2a ~ 3 { ~ c o c 2 + (a + l)c? } , 

^3 = c 3a - 4 {-^c 3 + (2 + 3a)coCiC2-(2 + 3a)(l + a)c?/2}, 
and so on. So, for 5 of (|3T29l) . x~ a = c^ l v{l + c 5(v a , x*)) &tv = l-u. Now apply (|330l) . □ 



Lemma 2.1 For in C k , 

k 

E\[{l-U n ^f =b n (r:G), (2.7) 

i=l 

where 

k 

b n (r : 0) = J] b (n - n-un - n + 1 : 0<) (2.8) 

i=i 

and 6(a, : 9 = B(a, (3 + 0) /B(a, (3) . Also in / TO) . 

fe 

B n (r) = JjB(ri-n_i,n-ri + l). (2.9) 

i=l 



Note 2.2 Since B(a, (3) = oo /or i?e/3 < ; /or |g. 7| ) to be finite we need n — r« + 1 + Red > 
/or 1 < i < k. 



Proof of Lemma 2.1 Set I k = LHS^XJ}) = J B n (u : r) [IiUC 1 ~ Ui) 9l dui ■ ■ ■ du k integrated 
over < m < ■ ■ ■ < u k < 1 by ([TTIj) . So, (|277jl . (123]) hold for fc = 1. Set Sj = (m - Ui_i)/(1 ~ 
iii_i). Then 

h= Cu^il-urt / 1 (n 2 -n i r- ri - 1 (l-n 2 r- r2 - 1+e2 o!n 2 / J B n (r), 

jO Jui 

which is the RHS (|2.7p with denominator replaced by the RHS (|2.8p . Putting = gives 
(12771) . (US} for = 2. Now use induction. □ 



Lemma 2.2 In Lemma 2.1, the restriction 

1 < r\ < • • • < rfc < n may 6e relaxed to 1 < r\ < ■ ■ ■ < r k < n. (2.10) 
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Proof For k = 2, the second factor in RHS (23]) is_6(r 2 - r x ,n - r 2 + I : 9 2 ) = f{9 2 )/f(0), 
where f{9 2 ) = T(n — r 2 + 1 + #2)/T(n — r\ + 1 + = 1 if r 2 = r i and the first factor is 
b(ri,n — r\ + 1 : 0\) = E(l — U nri ) ei . Similarly, if r-j = rj_i, the ith factor is 1 and the product 
of the others is E]J k j=1 .^(1 - J7 nrj )^, where 9* = 6j for j ^ i - 1 and 0j = 0j_i + (9; for 
j = i - 1. □ 

Corollary 2.1 In any formulas for Eg (X n ^ r ) for some function g, 112. 1 U\) holds. In particular 
it holds for the moments and cumulants of X nr . 

This result is very important as it means we can dispense with treating the 2 k ~ 1 cases 
(rj < r i+ i or r^ = r i+1 , 1 < i < k — 1 separately. For example, Hall (1978) treats the two cases 
for cos(X n)I ., X njS ) separately and David and Johnson (1954) treat the 2 fc_1 cases for the fcth 
order cumulants of X n ^ separately for k < 4. 

Theorem 2.2 Under the conditions of Theorem 2.1, 

k co 

EU X ti= E C n ^---C lk ^ k b n {v:\a-~6 /a) (2.11) 

i=l ii,...,i fe =0 

with b n as in \2. 8\) . All terms are finite if ReO < (s + l)a, where Sj = u — ri. 

Lemma 2.3 For a, (3 positive integers 9 in C , 

b(a,(3:9) = J] (1 + 6/j)~ 1 . (2.12) 

3=0 

So, for 6 in C k , 

k Sj_i 

b n (r:d) = H H (l + 6/j)- 1 , (2.13) 

1 = 1 j = Si + l 

where Si = n — r^ and ro = 0. 

Proof: LHS (I2TT21 = T(f3 + 9)T(a + I3)/{T((3 + 9 + a)T{(3)}. But F(a + x)/Y{x) = (x) a , so 
(I2T21) holds, and hence ([2TT3]) . □ 

From (|2,8p we have, interpreting Y\i=2 as 1 when k — 1, 

Lemma 2.4 For Si = n — r^, 

b n (r : 0) = B(s : 6)n\/T (n + 1 + h) , (2.14) 
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where 

k 

B(s : 6) = r (si + 1 + h) (siiy 1 J] fo - *i + 1 : #i) 

1=2 

does not depend on n for fixed s. 
Lemma 2.5 VFe have 

OO 

n!/r(n + 1 + 6) = rT e ^ ^)n~\ 

i=0 

where 

e o (0) = 1, ei(0) = -(0) 2 /2, e 2 (0) = (0) 3 (30 + l)/24, 

e 3 (0) = -(0) 4 (0) 2 /(4!2), e 4 (0) = (0) 5 (150 3 + 3O0 2 + 50 - 2)/(5!48), 

e 5 (0) = -(0) 6 (#) 2 (30 2 + 70 - 2)/(6!16), 

e 6 (0) = (0) 7 (630 5 + 3156) 4 + 3150 3 - 910 2 - 420 + 16)/(7!576), 
e 7 (9) = -(0) 8 (0) 2 (90 4 + 540 3 + 510 2 - 580 + 16)/(8!144). 

Proof: Apply equation (6.1.47) of Abramowitz and Stegun (1964) for i < 2 and Withers and 
Nadarajah (2007b) for i<7. □ 

So, (|2.1ip . (|2.14p yield the joint moments of X n ^n" l / a for fixed s as a power series in 
(l/n,n" a ): 

Corollary 2.2 We have 

k oo 

E [J X n,n- Si = E nT ( n + 1 + ^ " C i ( S : ^ ' ( 2 - 15 ) 

i=l j=0 

where if) = 0/ot and 

Cj ( s : ^) = E ' ' ' C ik^k B (s : ia - V>) : ii H h ifc = j'} . 

So, if s, are fixed as n — > oo and Re(0) < (s + l)a, 

oo 

LHS&IB = n fl E n "*" jae i O - C i ( s : VO • (2.16) 

ij=0 

If a is rational, say a = M/N then 

oo 

LHSQFEB = n* 1 n- m/N d m (s : -0) , (2.17) 

m=0 
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where 

d m (s:<4>) = ^2{e i (ja-^ 1 )C j (s:iP):iN + jM = m} 
= ^2 {e m -ja {ja - ipi) Cj (s : V) : O'j < m/a] 
if N = 1; so for d m to depend on c\ and not just Co we need m < M. 

Note 2.3 The following dimensional checks can be used throughout. By U.5\) . dimci = 
(dimX) a+l P . By \2. 6]) , dimC^ = (dimX) . Also dimxi = (dimX)~ a and dimd m (s : tp) = 
dimCj(s : ip) = (dimX) 01 . 

Note 2.4 The leading term in 12.16\) does not involve c\ so may be deduced from the multi- 
variate extreme value distribution that the law of X n ^ n _ Si , suitably normalized, tends to. The 
same is true of the leading terms of its cumulants. See Withers and Nadarajah (2007a) for 
details. 



The leading terms in (I2,16h are 

[{1 - rT x < fa > 2 /2} C (s : 0) + n - a C (s : V) + O (n" 2a °)] , 

where 

ao = min(a, 1), 
C (s : = c B(s : -■0), 

k 

C\(s : ip) = CQ 1_a_2 Ci ^2 rfjB (s : alj — V>) 

i=i 

and for Ij = i for i m = <5 m ,-, that is J, m = I(m < j). For k = 1, 

Cj(s:ip) = Cj^(s + l)ja-V) 

C (s : V) = 4(s + = cj/ < s 

Ci(s : V) = ^ -a-1 Cl(s + l) a _^, = Vc^" a_1 Ci/ < S >^- a ■ 

Set vr s (A) = b{s\ — S2,S2 + 1 : A) = rij=s2+i V(l + ^/i) f° r ^ an integer. For example, 
tt s (1) = (a 2 + l)/(si + 1) and tt s (-1) = si/a 2 . Then for A; = 2, 

C (S : Al) = C 2 X < Sl > 2A 1 7T S (-A) 

= Cq (si — S2 for A = 1 

= Cq < «2 - 2 >2 X < «2 >2 _1 for A = 2 

d(s : Al) = Xc 2 x - a - 1 c 1 < Sl >-l_ a {tt s (-A) + vr s (a - A)} 

= Ac^ a ci < si >2i a {si/s 2 + 7r s (a - 1)} for A = 1 
= \cl~ a ci < si >±\ {< si > 2 < s 2 >2 _1 +vr s (a - 2)} for A = 2. 



and 
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Set A = 1/a, Y ns = X ntn _ s /(nco) x and E c = Ac a 1 c±. Then for s > A — 1 

EY ns = {1 - n- 1 < A > 2 /2] < s >^ +n~ a E c < s >^ a +0 (n~ 2a °) (2.18) 

and for s% > 2A — 1, S2 > A — 1, s\ > s 2 , 

^y nSi y nS2 = {1 - n- 1 < 2A > 2 /2} £ 20 + n- a E c D a + (n- 2a °) , (2.19) 

where B 2 o =< si > 2A 1 ^(-A), D a =< si > 2 l_ a {^(-A) + vr s (a - A)} and 

Covar (Y nsl ,Y nS2 ) = F + Fi/n + £ c F 2 /n + (n~ 2a °) , (2.20) 

where F = B 20 - < si >^< s 2 > A ~ 1 , =< A > 2 < «i >^ X < s 2 >^ - < 2A > 2 £20/2 
and F2 = D a — < s\ >^ < s 2 >\- a ~ < s l > \- a < 52 >^ ■ Similarly, we may use (|2.16p 
to approximate higher order cumulants. If a = 1 this gives and Covar (Y nSl ,Y nS2 ) to 

0(n~ 2 ). 
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Example 2.1 Suppose a = 1. T/ien y ns = X nin _ s /(nco), E c = c a cq, B- 
(si - 1) _1 S2 S F =< si >2 1 s^ 1 , £> a =< si >^l a G a , where G a = sis^ 1 + ir s (a - 1) fi 
si > s 2 , G a = 2 for si = s 2 and F 2 = D a - s^ 1 < s 2 >J"i a -Sj 1 < s x >i\- So, 



or 



EY ns = s" 1 + n~ a E c < s > 1 -_ 1 a +O(n- 2ao ) (2.21) 
for s > and IKW\) - ^EM) hold if 

si > 1, s 2 > 0, si > s 2 . (2.22) 
^4 /itt/e calculation shows that Cq(s : 1) = CqB^q, C\(s : 1) = c^T a ~^c\B k ., and 

k 

E\\Y n , Sl = {l + n- 1 < k> 2 /2}B k0 + n- a E c B k + O(n- 2a0 ) 
i=i 

= m (s) + n~ 1 m 1 (s) + n~ a m a (s) + O(n~ 2ao ) 
say for Sj > k — i, 1 < i < k and s\ > • • • > s k , where 



B k . = ^ B k j, 

= IJl/Csi-fc + l), 



i=l 



-Bfcj = jQ (sj - A; + a + i) 1 < Sj - k + j + 1 > a _i |^[ (s, - k + i) 1 

i=l i=j+l 
k-1 

Bkk = Y\(si - k + a + i)~ l < s k >^ a 
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for Si > k — i and 1 < j < k. For example, B\q = s\, B20 = (si — 1) 1 s 2 and B 3 q = 
(si - 2) _1 (s 2 - 1) _1 S3 l . So, n n (s) = K(Y nsi ,. . . ,Y nSh ) is given by K n (s = k (s) + n~ x K;i(s) + 
n~ a K a (s)+0(n~ 2a °), where, for example, writing T, 3 a(si)b(s2S 3 ) = a(si)b(s2S 3 )+a(s2)b(s 3 si)+ 
a(s 3 )b(sis 2 ), 

3 3 
ko(sis 2 s 3 ) = m (sis 2 s 3 ) - y^m (si)mp (s 2 s 3 ) + 2]~[m (s^) 

8=1 

= 2(si + s 2 - 2)D(siS 2 S 3 ), 

3 

«i (sis 2 s 3 ) = mi (sis 2 s 3 ) - (si)mi (s 2 s 3 ) 

= 2 {s 2 (1 - 2si) + si - si } /D (sis 2 s 3 ) since mi(si) = 0, 
3 

K a (siS2S 3 ) = m a (sis 2 s 3 ) - ^ {m (ai) m a (s 2 s 3 ) + m a (si) m (s 2 s 3 )} 

3 

+2 ^ ra (si) m (s 2 ) m a (s 3 ) , 
w/iere -D(sis 2 s 3 ) =< s x > 3 < s 2 > 2 s 3 . 

Consider the case a = 1. TTtera K a (sis 2 s 3 ) = so 

k„ (sis 2 s 3 ) = 2 {si + s 2 - 2 + n- 1 (s 2 (1 - 2si) + si - sf ) } /£> (sis 2 s 3 ) + O (ra- 2 $2.23) 

Set s. = X^j=i s i- 

Bi. = B u - 1, B22 = l/s 2 , £ 22 = l/s 2 , B22 = si, 
B 2 . = + s 2 1 = ( s i + s s) / (sis 2 ) , 

£31 = (s 2 - S3 \ 5 32 = (si - sj 1 , B 33 = (si - l)" 1 S2 \ 

#3- = {S2 (S- - 2) - S 3 } (Si - < S 2 > 2 ' 1 S3 1 , 

B41 = (82 ~ 2)- 1 (S 3 - I)"' Sj 1 , 5 4 2 = (Si - 2)" 1 (S3 - I)" 1 Sj 1 , 

S 43 = (si - 2)" 1 (s 2 - I)" 1 sj\ S44 = (si - 2)- 1 (s 2 - I)" 1 S3 1 , 

B 4 . = {s.s 3 (s 2 - 2) + s 3 (s 2 - 4s 2 + 4) - s 2 s 4 } {(s x - 2) < s 2 - 2 > 2 < s 3 > 2 s 4 } _1 . 

^4/so E'c = Cq 2 ci, -D a = s^ 1 + «2 , i*b = 0, and 

EY ns = s~ l + n~ l E c + O (n~ 2 ) for s > 0, (2.24) 
EY njS1 Y UtS2 = (1 - n- 1 ) B 2 + n- l E c D a + O (n" 2 ) if (KM holds, (2.25) 
CWr (F njSl ,y niS2 ) = < s x >" 1 s^ 1 (s - tTV) + O (n^ 2 ) j/ B) /io/ds. (2.26) 

In the case a>2, \2.2J^ - \2~2b^i hold with E c replaced by 0. In the case a < 1, h2.19\) - [2~Jl\) 
with ao = a give terms 0(n~~ 2a ) with the n~ l terms disposable if a < 1/2. 

We now investigate what extra terms are needed to make (|2.24p - (|2.26p depend on c when 
a = 1 or 2. 
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Example 2.2 a = (5 = 1. Here, we fine the coefficients of n 2 . By {2.11 ), 

2 

d 2 (s :tp) = J2 e 2-j (j - $l) Cj(s :ip)+e 2 (-$i) C (s : V) 

+ei (l-^i)Ci(s:^) + C , 2 (s:^) 
= C 2 (s : V) if $x = l or 2. 

For k = 1, C 2 (s : i/O = C 2 ^(s + 1)2-0, w/iere C 2 ^ = ^cjj 4 {coc 2 + — 3)cf/2}, so d 2 (s : 
1) = (s + 1)F C; where F c = Cq 3 (coc 2 — cf), so in {2. 2$ we may replace 0(n~ 2 ) by n~ 2 (s + 
1)F cC q 1 + 0(n~ 3 ). For k = 2, 

C 2 (s:l) = ^{C il C j iB(s:0J-l):i + j = 2} 

= C 01 C 21 {B(s : 0, 1) + B(s : 0, -1)} + C 2 u B(s : 0), 

w/iere B(s : 0, A) = b(s\ - s 2 ,s 2 + 1 : A) = vr s (A), so d 2 (s : 1) = C 2 (s : 1) - D 2 ^ S H C + Cq 2 c^ ; 
where D 2s = (s 2 + l)(si + + sis^ 1 , H c = Cq 2 (cqc 2 — c 2 ) and in 112. 25\) we may replace 
0(n~ 2 ) by n~ 2 d 2 (s : l)cg 2 + 0(n~ 3 ). Upon simplifying this gives 

Covar (Y n)Sl ,Y njS2 ) =< si >^ s^ 1 (l - n _1 si) - Cq 2 H c F 3<s n~ 2 + O (n~ 2 ) , 

w/iere F 3jS = (s 2 + 1)/ < s\ > 2 +s 2 • 



Example 2.3 a = l, (3 = 2. So, a = 2, A = 1, ip = 6. By [BTfy , 

1 

3=0 

= e 2 (-^ 1 )C (s:^) + C 1 (s: V) 
= Ci(s : V) ififri = 0,1 or 2. 

For k = 1, 



so Fy ns = s" 1 + Cq 3 ci(s + l)n~ 2 + 0(n- 3 ) /or s > 0. For k = 2, Ci(s : 1) = c^ci-D^ /or 
L> 2s above, so 

EY ntSl Y n , S2 = (1 - n" 1 ) (si - I)" 1 s 2 1 + n" 2 C Q 3 Cl D 2 , s + O (n" 3 ) 

and 



Covar (Y n)Sl ,Y njS2 ) =< s\ > 2 1 s 2 1 (l - n 1 s 1 )-n 2 c 3 ciF 3jS + O 
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3 Examples 



Example 3.1 For Student's t distribution, X = t jy has density 

oo 

{l + x 2 /NY"g N = Y J d l x-^\ 



i=0 



w/iere 7 = (TV + l)/2, c/jv = r( 7 )/{ VA^r(iV/2)} and cfc = ^)N^ +i g N . So, tfjfl) holds with 
a = N, (3 = 2 and Ci = di/(N + 2i): 

c = iV 7 " 1 ^, 

ci = - 7 iV 7+1 (A r + 2)- 1 5J v = -iV 7+1 (iV + l)(iV + 2)- 1 5i v/2, 

c 2 = ( 7 ) 2 iV 7+2 (A r + 4)~ 1 5 7v/2, 

c 3 = -(7) 3 A r7+3 G JV (^ + 6)- 1 /6, 

and so on. So, a = 2/N and ( 2.17]) gives an expression in powers of n~ a l 2 if N is odd or n~ a 
if N is even. The first term in [2.17\j to involve c\, not just cq, is the coefficient of n~ a . 



Putting N = 1 we get 



Example 3.2 For the Cauchy distribution, hi. 5^ holds with a = 1, (3 = 2 and Cj = (— l) l (2i + 
l)" 1 ^- 1 . So, a = 2, V = 6, C Qil = = -^2-^/3, C 2il = ^ 4 "^{l/5 + - 5)/a} 

and C 3V , = -V>tt 6 -^{1/105 - 2^/15 + (ij> + 1) 2 /162}. By Example 2.3, Y ns = (■n/n)X n . n ^ s 
satisfies 

EY ns = s' 1 - n~V(s + 1) + O (n- 3 ) (3.27) 

for s > and when \2. 22\) holds 

EY n<sl Y n , S2 = (1 - n" 1 ) ( Sl - I)" 1 s^ 1 - n~VD 2iS /3 + O (n~ 3 ) (3.28) 

for D 2 , s = (s 2 + l)/(si + 1) + si/s 2 and 

Comr(y njSl ,K„ iS2 ) =< Sl >J X s^ 1 (l - n -1 ai) + n~ VF3 )S /3 + O (n^ 3 ) 

/or F 3iS = (s 2 + 1)/ < si >2 +s 2 l . Hall (1978, page 274) gave the first term in HJFfy 
and h3.28\) when s± = s 2 but his version of h3. 28\) for s\ > s 2 replaces (si — l) _1 s^ and 
D 2 , s by complicated expressions each with s± — s 2 terms. The joint order of order three for 
{Y n<Si , 1 < i < 3} is given by \2.23\) . Hall points out that F _1 (n) = cot(-7r — iru), so F~^(u) = 
IXot 1 ~ u) 2l ~ l C % i, where C a = {-47r 2 ) i 7r- 1 B 2 i/(2i)\. Note 2.1 could be used. We have not 
done so. 



Example 3.3 Consider the F distribution. For N,M > I, set v = M/N, 7 = (M + N)/2 
and g MN = v M / 2 /B(M/2, N/2). Then X = F M>N has density 

00 

x M l 2 (1 + vx)^g MN = v-<x- N ' 2 (1 + v ',- ') "g MN =Y J d i x- N,2 -\ 

i=0 
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where d* = /imjv(7V and h MN = qmnv' 1 = v~ N/2 / B{M/2, N/2). So, for N > 2, l2~6\) 
holds with a = N/2 — 1, (3 = 1 and Ci = di/(N/2 + i — 1). If N = 4 then a = 1 and Exam- 
ples 2.1-2.2 apply. Otherwise A2.18\) - [2.2(J\) give EY njS , EY n)S1 Y n ^ S2 and Covar(Y nsi ,Y n>S2 ) to 
O(n- 2ao ), where Y n , s = Y„ >n _ s /(rac )A, A = l/a, a = 2/(N - 2), a = min(a, 1) = a if N > 4 
and oq = min(a, 1) = 1 if N < 4. 



Example 3.4 Consider the stable laws. Feller (1966, page 549) proves that the general stable 
law of index a £ (0, 1) has density 

oo 

Ixl-^^OfcCa^), 

fc=i 

where afe(a, 7) = (l/7r)r(/ca + 1){(— sin{A:7r(7 — a)/2} and | 7 |< a. So, /or x > its 
distribution F satisfies 112. 6\) with (3 = a and q = aj + i(a, 7)7 _1 (i + 1) • Since a = 1 i/ie /irsi 
too moments ofY n ^ s = X„ in _ s /(nco) A , where A = l/a are gu>en to 0(n~ 2 ) \2.18\) - [2.20\) . 



Example 3.5 Finally, consider the second extreme value distribution. Suppose F{x) = exp 
(— x~ a ) for x > 0, where a > 0. Then U.5\) holds with (3 = a and C{ = (— l) 4 /(i + 1)!. Since 
a = 1 the first two moments ofY n)3 = X n , n - s /n 1 l a are given to 0(n~ 2 by ll2.18\) - [2. 20\) . 
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Appendix A: An Inversion Theorem 



Given xj = yj/jl for j > 1 set 

oo oo 

S = S(t,x) = Y,xjt j = S(t,y) =J2vjt j /jl (3.29) 
j"=i j"=i 

The partial ordinary and exponential Bell polynomials B ri (x) and B ri (y) are defined for r = 
0, 1, ... by 

oo oo 
r=i r=i 

So, B r o(x) = B r o(y) = 5 r o (1 or as r = or r / 0), B r i(\x) = \ % B r i(x) and B r i(Xy) = 
X-B^y). They are tabled on pages 307-309 of Comtet (1974) for r < 10 and 12. Note that 

oo oo 

{l + XS) a = Y / t r d r = Y, tr Cr/r\, (3.30) 

r=0 r=0 

where 

C r = C r (a, \,x) = J2 Bri{x) ( a \ X i (3.31) 
i=o \ l ' 

and 

r 

C r = C r (a, \,x) = ^2 B ri (y) < a >; A*. 
i=0 

So, Co = 1, Ci = aAxi, C2 = a\x2+ < a >2 X 2 x\/2, C3 = aXx^+ < a >2 A 2 xiX2+ < a >3 
A 3 a;f/6 and Co = 1, Cl = a\y\, C2 = a\y2+ < a >2 X 2 y\. Similarly, 

00 00 
log(l + XS) = t r D r = ^ t r D r / r \ 



r=l r=l 



and 



exp(AS) = 1 + t r 5 r = 1 + E * r -Br/r!, 



r=l r=l 



where 



5 r = D r (X,x) = -^2B ri (x)(-Xy/i\, 



D r = D r (X,y) = -J2B ri (y)(-Xy/(i - 1)!, 

i=i 
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B r = B r (X, x) = B ri {x)X i /i\ 

8=1 

and 



B r = B r {\y)=Y,B n {y)X l . 



i=l 

Here, B r (l,x) and B r (l,y) are known as the complete ordinary and exponential Bell polyno- 
mials. If Xj = yj = for j even, then S = t^ 1 Y^jLi Xjt 2 i , where Xj = x%j-\, so 

oo oo 

S* = t~* ^ t 2r B ri (X) and exp(AS) = 1 + ^ t k B k , 

r=i k=l 

where 

B k = {SkCXJAV*! : » = 2r - fe, jfc/2 < r < A;} . 
The following derives from Lagrange's inversion formula. 

Theorem A.l Let k be a positive integer and a any real number. Suppose 

oo oo 

v/ u = ^2x iU ia = ^2 y y a /i\ 

i=0 i=0 

with xq 7^ 0. Then 

oo oo 

(u/v) k = ^v ia = Y,y> ia /(ia)\, 
i=0 i=0 

where x* = x*(a, k, x) and y* = y*(a, k, y) are given by 

i 

x * = kn~ 1 d l (-n, 1/so, x) = fcxp n £(n + (-zo)~ j /j! (3.32) 

3=0 

and 

i 

y * = fcn-iCi (-n, l/y , y) = ^ £(n + l)i-i#y (y) (-yo)~ j , (3.33) 
respectively, where n = k + ai. 

Proof: u/u has a power series in v a so that (u/v) k does also. A little work shows that (|3.32|) - 
(|3.33p are correct for i = 0, 1, 2, 3 and so by induction that x*x 1 q and y*yg a are polynomials in 
a of degree i — 1. Hence, (|3.32[) - ()3.33j) will hold true for all a if they hold true for all positive 
integers a. Suppose then a is a positive integer. Since v/u = xq(1 + Xq 1 S') for S = S(u a , x) = 
S(u a ,y), the coefficient of u ai in (v/y)~ n is Xq n Ci{— n, 1/xq,x) = y^ ra Ci(— n, 1/yo, y)/(n — k)\. 
Now set n = k + ai and apply Theorem A in Comtet (1974, page 148) to v = f(u) = 

££o*iu 1+oi - □ 

Note A.l Comtet (1978, page 15, Theorem F) proves \3.32i) for the case k = 1 and a a 
positive integer. 
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